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In spin-lattice models with order parameter conserved, we generalize the idea of spin diffusion 
incorporating a variety of factors as possible driving forces, including external field and temperature. 
The Kawasaki dynamics in the Gaussian model and the one-dimensional Ising model are studied as 
specific examples. In the obtained diffusion equation, the term describing the diffusion induced by 
the inhomogeneity of the magnetization itself is unaffected and is believed to vanish near the critical 
point. Meanwhile the nonvanishing diffusion induced by the inhomogeneity of the environment may 
be coupled to the spin configuration and weakened by thermal noise. Interesting dynamic behavior 
is observed as a result of a competition of internal and external inhomogeneities and time scales. 
Several interesting examples are visualized, and the concept of local hysteresis is proposed in this 
spin-conserved dynamics. 
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I. INTRODUCTION 

In the past decade, great interest has been aroused by the response of a cooperative system to external perturbations, 
e.g. an external field varying with time or a scanning temperature (see Ref. [1-7] and references therein). In particular, 
significant effort has been made in the dynamic phase transition and hysteresis of the Ising model subject to an 
oscillating field [1], using Monte Carlo simulation and Mean-Field approximation based on Glauber's dynamics. Yet 
the dynamics with conservation law^, i.e. a model B in the classification of Hohenberg and Halperin [10], has been 
much less studied. Such dynamics may describe, for example, the phase separation of a binary mixture subject to 
inhomogcneous temperature or gravity. 

Generally speaking, in spin-lattice models with the order parameter conserved, it is possible to describe a diffusion 
process with a simple diffusion equation. In homogeneous environment, it is natural to guess that the spin flow is 
proportional to the gradient of the spin, and the equation will be 

'-^ = vvMr,t), (1) 

where q (r, t) denotes the local magnetization and T) is a diffusion coefficient. This equation has been exactly obtained 
in the Gaussian model and T) is found to vanish at the critical point (see Ref. [11], and the following section). It has 
also been studied in other systems, including the Ising model [12]. 

Diffusion is always a result of inhomogeneities. The diffusion process characterized by Eq. (1) can be viewed as 
being caused by the inhomogeneity of the magnetization itself, while diffusion driven by external inhomogeneities 
remains to be an unsolved interesting problem. Many studies on the dynamic response of the Ising model deal with 
a uniform field varying with time [1]. However, such a field has no effect on the Kawasaki dynamics [12] with spin 
conservation. So here is a more difficult task, which should incorporate both spatial and temporal inhomogeneities. 
In relevant studies several interesting problems have been treated, e.g., the random-field spin-conserving Ising model 
[13-16] and, more recently, the domain growth in a temporally constant field varying linearly in space [17,18]. In the 
present study, we narrow our scope to the Gaussian model and the one-dimensional Ising model, and focus on the 
modification of the above diffusion equation by the inhomogeneity of field and temperature. 

As is well known, the Ising model has wide and important applications in various fields. Unfortunately, as we shall 
' see below, the obtained diffusion equation can not be solved. Meanwhile, it is analytically tractable in an idealization 
of the Ising model, the kinetic Gaussian model, which may serve as a starting point and a test field of new ideas. 
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'^In the absence of driving force, the Ising model with spin conserved is relatively well understood: Quenched from a high 
temperature disordered state, some kind of order rises with the appearance of coarsening and competing up-and-down domains, 
of which the length scale grows as a power law with time [8,9]. 
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In this article, the rigorous treatment of the Gaussian model will be combined with the less quantitative discussions 
and Monte Carlo simulations of the Ising model, along the same line, for a comparison. Section II deals with the 
diffusion induced by inhomogeneous external field, and in Sec. Ill we extend the scope to inhomogeneous temperature. 
(Although here we treat these two problems separately, it is possible to put them together.) Section IV is the 
summarization with some discussions. 



II. INHOMOGENEITY OF THE EXTERNAL FIELD 



In this section, based on the master equation, wc shall study diffusion process in an inhomogeneous external field. 
Such a process is subject to two factors acting together, the inhomogeneity of the field and that of the magnetization 
itself. This competition may lead to interesting and rich dynamic behavior, as we shall sec below. Before wc move to 
the specific models, we briefly review the governing dynamic mechanism, spin-pair redistribution mechanism, which 
is a natural generalization of Kawasaki's spin-pair exchange mechanism [12]. 

In the Ising model, spins can only take two opposite values, +1 and —1, and spin-exchange is simply a natural 
choice; but it turns to be not as capable when applied to other more complex systems. A natural generalization is as 
follows: a pair of nearest neighbors selected at random, at and aj, may take any values, ai and aj, as long as their sum 
keeps conserved, and this is called spin-pair redistribution (clearly in the Ising model it reduces to the spin-exchange). 
We list several important equations in Appendix A, the details of which can be found in Ref. [11]. 



A. The Gaussian model 



First we will study the three-dimensional kinetic Gaussian model in an inhomogeneous external field, which might 
also be varying with time (e.g., a traveling wave). In lattice models, a traditional way is to assign to each vertex, cTj, 
a reduced field, Hi (t) [1]. In the kinetic Gaussian model, 

-/3H = if^aia, + ^^i/,(t)a„ (2) 

where K = J/ksT, and J is the coupling strength. The spins take continuous value from — oo to -|-oo, and the 
probability of finding a given spin between c7j to Gi+dcji is assumed to be the Gaussian- type distribution / (a-,) dcji = 

•y'^^exp (— 517^) daj . where h is a distribution constant independent of temperature. Thus, a summation over the 

spin values will be in fact an integration, — !■ / (c) da. 

Now, we substitute the Hamiltonian (2) to the evolving equation of single spin, Eq. (A3). A linear equation can 
be obtained (three-dimensional) 

+ {qij+i,k — 2qijk + q-Lj-i,k) + {qi,j,k+i — '^q-ijk + q-Lj,k-i)] 

+K [2 {2qi-ij^k — qi-i,j+i,k — qi-i,j~i,k) + (2gi-i,j,fe — qijk — qi-2,j,k) 
+2 (2gj+ij_fc — q'i+ij+i^fe — j_i_fe) -I- {2qi+ij^k — qijk — qi+2,j,k) 

+2 {^qi,j-l,k — qi,j-l,k+l — 9i,j-l,fe-l) + {^qi,j-l,k — qijk — qi,j-2,k) 

+2 (2gi,j+i,fe — qi,j+i,k+i — qi,j+i,k-i) + (2gij+i,fe — qijk — qi,j+2,k) 

+2 {2qi.j^k-i — qi+i,j,k-i — qi-i,j,k-i) + (2gi,j,fe-i — qijk — qi,j,k-2) 
+2 (2(7ij\fe+l — qi+ij^k+l — qi-l,j,k+l) + {'^qi,j,k+i — qijk — qi,j,k+2)] 

1 

[(2-ffijfc - Hi+ij^k - j,fe) 

+ i^Hijk — Hij+i^k — -ffi,j-i,fe) + {'^Hijh — Hij^k+1 — Hij^k-i)]} (3) 

This equation is in fact not as complex as it seems to be. Each term in a bracket is a second-order derivative either of 
the magnetization or of the external field, and they will cancel each other if a summation is taken, guaranteeing the 
conservation of spin. With lattice constant a we can transform the above equation (and similar equations for lower 
dimensions) into 
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(4) 



where D is the dimensionahty. It reveals an important feature that, surprisingly, the field and the spins are not coupled, 
which is mainly a result of the integration taken from — oo to +00. The influence of the inhomogeneous external fleld 
rigorously takes the form of a second order derivative, with a prefactor getting weaker at higher temperature (weakened 
by thermal fluctuations). 

Discussions: (1) If the external field is spatially homogeneous, the system behavior can be described by a simple 

diffusion equation [11]. 

(2) By contrast, if the external field varies solely in space, an equilibrium state can be obtained by setting 
dq{r,t)/dt = 0. Wc find that V^g (r) oc V'^H{r), except at the critical point Kc = b/2D. The susceptibility, 
X ~ V^g (r) /V^i? (r), can be obtained. To visualize this process, we suppose that the shape of the external field is 
a Gaussian packet: Above the critical point, the magnetization will reach the stable equilibrium also in the shape 
of a Gaussian packet, and x is finite and positive. As the system is cooling and approaching the critical point, we 
shall observe that the peak of the magnetization (as well as x) tends to positive infinity. With the scaling hypothesis 
X ~ \T — Tc\~^ , the critical exponent 7 = 1 can be (in fact generally) obtained. When the temperature is below 
Tc, the equilibrium value of the magnetization (and x) will suddenly become negative infinity. As the temperature 
continues to decrease, the magnetization becomes flatter, and x is finite and negative. In this region the equilibrium 
is obviously unstable. 

A comparison of the two dynamic versions, Kawasaki's and Glauber's, of the Gaussian model shall be interesting. 
In Glauber's version, we can obtain [20] 

|(<7,(i)) = -{am + IY1 Kij{aj{t)) + (t) (5) 

and 

where M = J2k (^fe) /-^' ^^'^ ^ — Sfe Hk/N. The external field does not change the critical point either, and we can 
also get the equilibrium magnetization by setting dM/dt = 0. The situation here is quite similar: Above the critical 
point, the equilibrium is stable, and it becomes unstable below the critical point, and x changes sign. In other words, 
in both versions, the observable state below is a fast growing dynamic process. 

(3) Prom Eq. (4) we are able to obtain the whole dynamic process, and here we discuss two interesting one- 
dimensional examples. 

The first example: We assume a time-independent field H = Hq cos (kx) , and the initial magnetization is zero. 
Substitution into Eq. (4) yields a solution, 

q {x, t) = — {1- exp {-Vk'^t}) Hq cos {kx) , (7) 

where V = {b/2 - K) / {b + K), and A = aV [2 {b + K) ksT]. When K < = 6/2, V is positive, and q wfll 
approach the equilibrium value, q^q = {A/V) Hocos{kx), with the relaxation time r = l/Vk"^. A typical critical 
slowing down phenomenon is obvious when K — » K~ , since r — > 00. Interestingly, the speed, dq/dt ~ exp {—Vk'^t), 
does not really "slow down" because of the diverging prefactor A/T), and this is against our common knowledge about 
this phenomenon. 

The second example: We assume a time-dependent field H = HQCOs{kx — cot), which is a traveling wave, and we 
obtain a solution of Eq. (4) as 

Ak'^ 

q {x, t) = —^====Ho cos {kx — ojt + (p) , ip = arccot {Vk^/co) . (8) 

It describes a spin wave, which lags in phase with a temperature-dependent factor p. As a counterpart of the hysteresis 
loop studied previously, formed by the average magnetization and the spatially uniform field, here we introduce the 
concept of local hysteresis loop, which is formed by local magnetization and local field. It is a useful tool characterizing 
the response to oscillating perturbations in a simple way. Generally, the local hysteresis loop may vary in space, while 
in this specific example it has an elliptical shape that is spatially uniform and varies with temperature in the way 
illustrated in Fig. 1. 
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B. The Ising model 



Now we will turn to study the Ising model. The Hamiltonian is of the same form as Eq. (2). In the one-dimensional 
version of Eq. (A3), we first treat the two combined terms 



(J/jC ^ -f- (jfcilt; ^ 



i (cTfc + CTfcil) , li (Tk = (Tfc±l 

(1 - C7fecrfe±i)tanh if (crfe=pi - (Tfc±2) + (-fffc - Hk±i) 



, IftTfc 



We assume that the variance of the external field between two nearest-neighboring vertices is small (i.e. when the 
wavelength is much longer than the lattice constant), then 



tanh 



^ (crfc:^ 1 - CTfc±2) tanh2ii' + {Hk - Hk±i) 



!--(!- (Jk^iak±2) tanh 2K 



Substituting them into Eq. (A3), we get the evolving equation of single spins. 



d 
di 



ft = y [ilk+i -Qk) /a - {qk - Qk-i) /a] /a 

— ^{[(9fc+2 -%-i)/3a- {qk+i - ^-2) /3a] /a}tanh2ii' 



{[{{'^kO-k+io-k+2) - {ak-io-ko-k+i}) / a - {{ak-iako-k+i) - {akO-k-ia-k-2)) /a] /a}taxih.2K 



1 1 



2 fcsT 
1 1 



2kBT 



{Hk+i - Hk) 
{Hk - Hk-i) 



1 - {f^k^^k+i} - 2 ((^ ~ o-feO-fc+i) (1 - o-fe-iO-fe+2)) tanh^ 2K 
1 - {(TkCTk-i) - ^ ((1 - o-fe(7fe_i) (1 - <Jk+i(Jk-2)) tanh^ 2K 



If the external field is spatially uniform, we use (cfc'^^^ to denote {(Tk-i<Jk<^k+i) ■, and 



dqk 



M 



In Ref. [12], Kawasaki has used local equilibrium approximation and concluded that 



V 



M 



(9) 



(10) 



where T) oc 1/x is the diffusion coefficient vanishing near the critical point. If this is correct in the Ising model 
(tanh2ii'c = 1 for the one- dimensional case), it requires. 



fq 



(11) 



to be true at least near the critical point, otherwise the diffusion process in the Ising model cannot be described by 
such a simple equation. 

In the following wc will incorporate the inhomogeneity of the external field in our discussion. Let the function 
f {ka±a/2) denote the whole term 



/ {ka ± a/2) = 1 - {<Jk(Jk±i) - 2 ((1 ~ ^k(^k±i) (1 - crkT-i-^k±2)) tanh 2K 
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and we get 



dq fdq\ 1 a2 f^,,d^H OH df{x) 



m \dt J 2A;sTV 

Compared with the Gaussian model, we find both similarities and difference. First, the diffusion term induced by 
the inhomogeneity of the magnetization itself retains the same form as Eq. (9), and this is also true for the Gaussian 
model. Second, the influence of the external field is also weakened by a prefactor 1/T. However, here the field is 
coupled to the spins. And, a somewhat surprising prediction of Eq. (12) is that the driving force may partly come 
from the gradient of the field, which is at the same time coupled to the inhomogeneity of the spin configuration. For 
example, there will still be field-induced "diffusion" if the field varies with space linearly and thus d^H/dx^ — 0, 
provided that the magnetization is inhomogeneous and thus df /dx ^ 0, and this case has been particularly studied 
in Ref. [17,18]. Beside this case, a lot of issues will make us even more curious, including, as studied in the Gaussian 
model, the effect of a travelling wave: Shall there also be a similar spin wave with the same time period and a phase 
lag? Shall there be new dynamic phases? What will be the shape of the local hysteresis loop? Here we provide a 
tentative answer with Monte Carlo (MC) simulations. 

In a one-dimensional model of N vertices with periodic boundary condition, different initial conditions are chosen 
while the total spin is kept around zero. We randomly pick a pair of nearest neighboring spins, Uj and Uj, and 
exchange them with a probability 

Wij = min {1, exp {-AEij/kBT)} , 

where AEij is the change in energy if ai and aj are exchanged. A MC step consists of N such actions. Time is 
measured in MC steps and the system is assumed to be of unity length. A travelling wave, H = Hq cos {kx — uit) is 
applied. 

Fig. 2 shows the local hysteresis loops of the 1st spin in a 100-spin chain at several typical temperatures, with 
oj = 27r/1000, k — 2tt and Hq = 4 (in units of J = ks)- (The shape of a loop is found to be independent of the position 
of the spin it describes.) Each point on a loop actually denotes the expected value of the spin, i.e., — p-, where 
is the probability that the spin takes upward direction and ]7_ the contrary. The results are to a great degree similar 
to the rigorous ones obtained in the Gaussian model. At K = 0.001 (Fig. 2(a)), the spin fluctuates around zero 
(actually it may be a very narrow ellipse with a tiny phase lag). At if = 0.05 (Fig. 2(b)), the loop is a tilted ellipse 
with a finite phase lag between and 7r/2. At about K = 0.4 ~ 0.6 {K = 0.5 in Fig. 2(c)), the loop evolves into a 
state that has two symmetries. It is symmetric about the horizontal axis, and the average magnetization over a period 
is zero. At the same time it is also symmetric about the vertical axis, and this indicates a phase lag of approximately 
7r/2. In the Gaussian model it agrees with the critical point, and here it also corresponds to a region where a dynamic 
symmetry loss occurs. First, the loss of the vertical symmetry: as temperature continues to decrease, the loops begin 
to transform into a shape which we are familiar with {K = 2.5 in Fig. 2(d)). It means that, though the spin value 
still varies with the same period, it no longer varies sinusoidally. Second, the loss of the horizontal symmetry: for a 
given spin, the average value over a period may increasingly deviate from zero as temperature decreases, while the 
extent and direction of the deviation are influenced by fluctuations and initial conditions. This may be accepted as 
an inherent property of the one-dimensional Ising model and is observed in the quasi-static limit. We find that the 
temperature at which a noticeable symmetry loss occurs is lowered by the magnetic wave. 

Here, we tentatively reveal the interesting features that result from a competition of the external and the internal 
inhomogeneities and time scales. A thorough study is needed to answer questions such as: What is the effect 
of wave speed, frequency and amplitude? What is the relationship between the area of the local hysteresis loop 
and the parameters (previous works in non-conserved dynamics have revealed several interesting scaling laws)? In 
higher dimensions, the effect of a magnetic wave on the formation and destruction of domains also deserves further 
investigation. In the case of a plane wave, we may expect to observe two coexisting length scales, one parallel to the 
wave direction, and the other perpendicular to it. This requires extensive MC simulations which are beyond the scope 
of this article. 



III. INHOMOGENEITY OF THE TEMPERATURE 



In Sec. II, we have investigated the field-induced diffusion. As mentioncxl in the Introduction, the inhomogeneity, as 
the driving force of the diffusion, may be actually of a broad meaning. In this section we shall focus on the influence of 
inhomogeneous temperature. Inspired by the dynamic response to a periodically altered temperature experimentally 
observed in ferroelectric systems [4] and charge-density- wave systems [5], there has been theoretical effort in simple 
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spin models [6,7]. With spin-nonconscrvcd dynamics, interesting behavior and thermal hysteresis have been reported. 
In the following we present our treatment of the Kawasaki dynamics of the Gaussian model and the Ising model. 

For simplicity we limit our scope to an one-dimensional model that consists of N spins. There is no (or homogeneous) 
external field and each vertex, ct^, is in contact with a heat reservoir of temperature T^. Thus the system's effective 
Hamiltonian 

-Pn^ff ({a}) = J2 ^^^^i^i- (13) 



where Ki = J/ksTi. With 



and 



(iJ> 



VK VT 



K T ("> 



we shall treat K for convenience in the following. 



A. The Gaussian model 

In the Gaussian model, we substitute the new Hamiltonian, Eq. (13), into the evolving equation of single spin, Eq. 

(A3), and get 



j^qk [t) -- qk+ 2[b+{Kk + Kk+i)/2] 



Substituting 



Kk + Kk-i 



Qk-i 



^ I (Qk + Qk+i) - I ^ ) qk+2 + b{qk+ qk+i) 



+ 



1 



2[b+{Kk + Kk-i)/2] 
Kk + Kk-i 



Kk + Kk+i 



qk+i 



{qk + qk-i) 



Kk-i + Kk-2 



qk-2 + b {qk + qk-i) 



(16) 



and 



into Eq. (16), we get 



Kk±i 



dK 



,d^K 



Kk ± + ^ci^^, Kk±2 ~ iffc ± 2a— ^ +2a^ 

ax'' ax ax'' 



dK 
dx 



qk±i « gfe ± a 



dq 
dx 



1 2 d.'^q dq n ^"^9 

-^a qk±2 ~ qk± 2a— + 2a 



dx^' 



ldq_ 2 36/2 -g^ /I d^q 
qdi ~ li' yqlk^ 



a2(2/s:' + 36) [ldq\ ( 1 dK 



+a' 



1 dK' 
'K'~d^ 



2K 

1 d^K' 
It dx^ 



q dx 



K' dx 



(17) 



(18) 



(19) 



where K' = K + b. The first term on the right hand side denotes the diffusion induced by the inhomogeneity of the 
magnetization itself, which we are familiar with. The influence of the inhomogeneity of K is described by the other 
three terms. The temperatiirc and spin are partially coupled. It is well established that the magnetization-induced 
diffusion will vanish near the critical point (in the one-dimensional model K^ = b/2), but the K-induced diffusion 
will not {neither will the field-induced one). With Eqs. (14) and (15), we can see that the role played by the thermal 
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noise is more difficult to analyze in i^-induced diffusion, compared with the case of field-induced diffusion. Besides, 
the third term on the right hand side is always positive, yet, surprisingly, this is still not against the conservation 
of the order-parameter. We can use a simple example to show it. Assume that at one moment the magnetization is 
homogeneous, and we have 



dq , 2 



1 dK'\^ ( 1 ^K' 



K' dx \K' dx^ 



The previous studies cited above mainly considered non-conserved dynamics when the temperature is varying with 
time but spatially homogeneous. Although a similar study of local thermal hysteresis is possible, here we turn to 
treat the contrary situation: spatially modulated but temporally fixed temperature. In the following we visualize the 
equilibrium state (if there is one) and the evolving process in three typical examples. The numerical simulations are 
based on Eq. (16), and the system consists of 100 spins with periodic boundary condition. We set 6 = 1 and assume 
that when t = 0, all the spins take the value of unity. 

The first example: The overall system is above the critical point {K < Kc ~ 1/2) and Ki = 1/4+ (1/8) sin (2i7r/100), 
i = 1,2, ...,100. The evolving process is shown in Fig. 3. We can clearly see that the system is approaching 
an equilibrium state and is slowing down as time passes. Changing the parameters, but keeping qi (0) = 1 and 
-K^max < Kc, we find that, interestingly, the spins keep to be positive in the evolution. 

The second example: The whole system is below the critical point and Ki = 1 + (1/8) sin (2i7r/100), i = 1,2, 100. 
Here we may expect the system to be fast growing. The results show an interesting self-organizing process, which can 
be characterized by two succeeding dynamic phases, a steady one, and a growing one. Fig. 4(a) shows a typical spin 
configuration during the period between t = and approximately t = 400. The magnetization takes a sinusoidal-like 
shape and the amplitude is growing very slowly. After the amplitude reaches a threshold value, we can observe the 
self-organization as shown in Fig. 4(b). Note that part of the magnetization, to gioo, still seems to be smooth. In 
fact, self-organization also occurs in this part, only later and weaker. After this period, the system moves into a fast 
growing phase with a fixed shape, as shown in Fig. 4(c). The contrast in the intensity of motion exhibited by the two 
regions, qi to (749 and q^o to gioo, persists if we change the temperature to, for example, Ki = 2 + (1/8) sin (2i7r/100). 
So it may be caused not by the absolute value of the temperature, but by the mutual influence between the relatively 
high temperature region and the relatively low temperature one. 

The third example: Here we study a rather interesting case, where Ki = 1/2 + (1/4) sin (2i7r/100). We can divide 
the system into two regions, qi g'50 and 951 gioo- In the first region, the observed evolution is similar to the 
second example, but with less peaks and lower growing speed (it is because the temperature is higher here). In the 
second region, however, the system behavior is similar to the first example, i.e. the magnetization is smooth and 
approaching equilibrium. We do not observe any strange phenomenon at K50 = Kc = 1/2. This is not surprising, 
and can be predicted with Eq. (19). 

B. The Ising model 

Now we turn to treat the one-dimensional Ising model. Similarly, by substituting Eq. (13) into Eq. (A3), we get 

1, ^ ,^1,. ^ f K^,+Kk K±i + Kk±2 \ 
= 2 V'^fc + ^fc±i) + 2 ~ ^kO'k±i) tanh I ak^i crk±2 1 

1 1 

= 2 ('^'^ + ^'^ii) + 4 ~ CTfcO-fcii) {akzfi - o-fe±2) tanh [(isTzpi +Kk + K±i + Kk±2) /2] 

T J (1 - cTfe(7fe±i) (cTfe^i + (7fe±2) tanh {{K^t + Kk - K±i - Kk±2) /2] . (20) 

Applying Eqs. (17)-(18) and similar equations for a^^\ we obtain 

dq fdq\ ^ 2ri , i,2,„j^,i 5X /15a(3) 3 dq 



odK fda^^^ dq\ ^d^K 
+« 7^ +0 



dx \ dx dx J dx^ 



(a(^) - g) . (21) 
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Compared with the rigorous results of the Gaussian model, there are similarities and differences. Similar to the 
Gaussian model, the first term at the right hand side, {dq/dt)j^j, is the diffusion induced by the magnetization itself 
(see Eq. (9)). It will vanish near the zero-temperature critical point. Besides that we can observe the nonvanishing 
iiT-diffusion, which will, in contrast, get infinitely sensitive near this point. A major difference is that, the terms of 
Eq. (21) do not take such relative forms as in Eq. (19). Thus, although the influence of dK/dx and d^K/dx^ still 
depend on the spin configuration, the dependence is in a different way. The system will get less sensitive when the 
temperature is increasing, and we may attribute this to the increasing infiuence of the thermal noise. 

IV. SUMMARY 

The present research is based on the Kawasaki-type spin-pair redistribution mechanism. We generalize the idea of 
spin diffusion incorporating a variety of factors as possible driving forces, including external field and temperature. 
Two models are selected for a detailed investigation: the Ising model, which is widely applicable, and the Gaussian 
model, which is idealized but provides a good basis for an analytical treatment. It is found that the two models are 
similar in principle, and the features are likely to be shared by other models and experimental systems. 

Generally speaking, the difl[usion equation can be written as 

dq f dq\ ( dq 



The right-hand side of the equation consists of two juxtaposed terms: One describes the diffusion induced by the 

inhomogeneity of the magnetization itself. It retains the same form as that obtained without any external inhomo- 
geneities, and is believed to be vanishing near the critical point. The other one is the "general diffusion" induced by 
the inhomogeneity of the environment, which may be coupled to the spin and contain both first and second derivatives. 
The latter, of which a more appropriate name is environment- induced self-organization, generally does not vanish near 
the critical point, and strongly depends on spin configuration and may be weakened by thermal noise. The concept 
of local hysteresis is proposed in this spin-conserved dynamics as a convenient tool to characterize the response of the 
system to oscillating external perturbations. The response of both the Gaussian model and the one-dimensional Ising 
model to a travelling electromagnetic wave is specifically studied. 

About the Kawasaki dynamics with external inhomogeneities, there remain a number of questions worthy of further 
investigations. For example, the domain growth in an Ising model subject to a travelling electro-magnetic wave^, the 
response to the temperature varying with time, and the dynamics in other models with continuous symmetry, such 
as the XY model [3] and the Heisenberg model, etc. We hope the interesting features revealed in this article may 
stimulate future research, probably with extensive numerical simulations. 
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APPENDIX A: THE SPIN-PAIR REDISTRIBUTION MECHANISM 

In spin-pair redistribution mechanism [11,19], two neighboring spins, cTjcr;, may change to any possible values, aj&i, 
as long as their sum are conserved. The master equation is 

|p({a}, i) = E E [-^^-^ ^^o<^i ^ ^^"^o p ({^} ; *) 

{jl) 

+Wji {aj&i ajai) P {{aijtj , aijtk} , aj,ai;t)] . (Al) 
The redistribution probability is in a normalized form determined by a heat Boltzmann factor. 



^Such processes are subject to several factors acting together: the inherent properties of the model, the moduling effect of the 
field, and a competion of time scales, namely the relaxation time of the system and the time period of the local field. 
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Wji {ajai aj&i) = —Saj+a,,aj+a, exp 



where the normahzation factor Qji is 



"3 



For single spins, the time expectation, g/j = (crfc) = '^fe-P({'''} i^)) is 

d 



dt 



qk{t) = -2Dqk W + 

w {a} 



(A2) 



(A3) 



where £> is the spatial dimensionality and denotes a summation taken over the nearest neighbors. 
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Caption of figures 



Fig. 1. Schematic plots of typical local hysteresis loops in the Gaussian model. 

Fig. 2. Typical local hysteresis loops for different temperatures in the Ising model. System size A'' = 100, and 
parameters w = 27r/1000, k = 2n and Hq = 4. 

Fig. 3. The evolution of the system above the critical point. The magnetization is approaching the equilibrium 
with its shape similar to the external field. Squares correspond to the system being at f = 50s, circles correspond to 
t = 300.S, and triangles correspond to t = 1000s. 

Fig. 4. The evolution of the system below the critical point, (a) The magnetization at t = 350s is in the steady 
phase and keeps a sinusoidal shape while growing very slowly, (b) The self-organization begins, and the magnetization 
at f = 400s is shown, (c) The system is in the growing phase and the magnetization is fast growing while keeping an 
fixed shape. Squares correspond tot = 500s, circles correspond iot = 510s, and triangles correspond tot = 520s. 
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